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We illustrate how Hawking's radiance from a Schwarzschild black hole 
is modified by the electrostatic self-interaction of the emitted charged parti- 
cles. A W.K.B approximation shows that the probability for a self-interacting 
charged particle to propagate from the interior to the exterior of the horizon is 
increased relative to the corresponding probability for neutral particles. We 
also demonstrate how the electric potential of a charged test object in the 
black hole's vicinity gives rise to pair creation. We analyze this phenomenon 
semiclassically by considering the existence of the appropriate Klein region. 
Finally we discuss the possible energy source for the process. 



I. INTRODUCTION 

The aim of this work is to study various effects that may take place when a charged 
object propagates on a background Schwarzschild spacetime. As elucidated by Vilenkin |J 
and corroborated by Smith and Will |2j and Lohiya ||, in contrast to the situation in flat 
spacetime, in the presence of a Schwarzschild black hole the self-field of the object, makes 
a nontrivial contribution to the object's energy as measured at infinity. Alternatively, this 
phenomenon may be interpreted as black hole polarization; the test object polarize the black 
hole so that image charges are formed on the black hole horizon. It was shown by various 
authors M that this effect is of unquestionable importance in black hole thermodynamics. 
One may think of two additional distinct effects that may take place on account of black hole 
polarization. The first is the emission of self-interacting charged particles. Here the black 
hole is isolated. The temperature of the black hole is assumed to be high enough so that 
the lightest massive charged particles may be emitted. The second effect is pair creation 
due to vacuum polarization. The black hole is allowed to interact with an exterior charged 
test object, and the potential generated by the charge is taken into account. Both processes 
seems to have the effect of increasing the luminosity of the black hole in a sense that will be 
clarified below. 

In the following sections we investigate these two effects. First we show that once the 
temperature of the Schwarzschild black hole ~ 1/M, is equal or greater than the rest-mass 
fi of a self- interacting charged particle ( fi < 1/M), the emission of these particles is more 
probable than the emission of similar neutral particles by a factor exp(7re 2 ), where e is the 
charge (unless the contrary is stated we use "natural units" in which h — c — G — k — 1 
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throughout). This is a direct result of the self-interaction of the charged particles. In analogy 
with the increased thermionic emission from a hot filament placed in a repelling electrostatic 
field, which has the effect of reducing the potential barrier at the edge of the metal, the 
effective electrostatic self-field experienced by the charges has the effect of increasing the 
emission ||. In order to demonstrate this we compute in Sec.[TI]the probabilities for finding 
particles in the exterior and interior of the horizon using the appropriate wave equation 
coupled to an effective electromagnetic self-field. It turns out that for a Schwarzschild black 
hole, the chemical potential is lowered (for both signs of charge ) by e 2 /8M. As pointed 
out before by several researches this suggests that the emission of charged particles 

would increase relative to the emission of neutral particles. We calculate an approximation 
to the absorption coefficient ("greybody factor") for the problem. We find that when the 
energy of the emitted particles is close to e 2 /8M, the absorption coefficient is of the form 
a(u, /i, e 2 /8M)(u; — e 2 /8M) plus small corrections. This means that in the energy range 
H < ijj < e 2 /8M, the black hole should start to superradiate . 

Thus far we have discussed the emission of charged self- interacting particles. In Sec.|II 
we change the settings and discuss the consequences of placing a charged test object in the 
exterior of a Schwarzschild black-hole. We illustrate how the field generated by the charged 
particle gives rise to pair production. We analyze the problem semiclassically and find the 
energy range for the existence of a Klein region. It was previously thought that this energy 
range does not exist in the case of an isolated nonrotating neutral black-hole || . 

In Sec. IV we address the problem of possible energy sources for the process. Ther- 
modynamical arguments suggests that if the black-hole is too massive, the entropy carried 
away by the pairs is not sufficient to compensate for the entropy lost by the black-hole. 
Accordingly we look into the dynamics of the test object, and investigate the role it may 
play in the supply of energy for the process. Finally in Sec. V we summaries our finidng. 



II. PARTICLE EMISSION RATE. 



We take now a semiclassical approach to the calculation of the transmissions probabili- 
ties. We will use the complex path method which was recently advanced by various authors 
|7^f§. Specifically Srinivasan and Padmanabhan |7j proposed a derivation of Hawking ra- 
diation without using the Kruskal extension as in |I0| . The point, of course, is that the 
Schwarzschild coordinates possess a coordinate singularity at the horizon. This bad behav- 
ior of the coordinates manifests itself as singularity in the expression for the semiclassical 
propagator near the horizon and a specific prescription to bypass it must be provided. The 
action functional may be constructed using the Hamilton- Jacobi method in the appropriate 
coordinates. It will be shown how the analytic continuation used in complex path method 
gives the result that the probability for particles to be found in the exterior of the horizon 
is not the same as the probability for particles to be found in the interior of the horizon. 
The ratio between these probabilities is of the form, -Pmterior = X ex P( — P £ ) -^exterior, where 
e is the energy of the particles and j3 = 8nM is the inverse Hawking temperature. \ is a 
factor of the form \ = exp(7re 2 ), where e is the charge of the test particle. This implies that 
the usual thermal distribution of neutral particles is modified by the addition of what may 
be interpreted as a charge dependent chemical potential. 
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A. Transmission probabilities 



In a Schwarzschild spacetime consider a minimally coupled test scalar field $ with mass 
[A coupled to the electromagnetic self-field A a generated by the charge e. The scalar field 
propagates in the metric 



ds 2 = -(1 - 2M/r)dt 2 + (1 - 2M/r)- 1 dr 2 + r 2 dfl 2 , 
where dfl 2 = d9 2 + sin^ 2 ^ 2 , and satisfies 



V, v - -eA a ) { V 



A a , linear in e, is the self-potential whose source is the object itself |l]-f|. 

eM , 



A,, 



2r 2 «■ 



(2.1) 



(2.2) 



(2.3) 



A a is divergence free accurate to 0(e 2 ). Using this and expanding the left hand side of 
Eq. (|2.2|), we obtain, 



\d r [r 2 (l - 2M/r)9 r $] - (1 - 2M/r)~ 1 (d t - % -eA^j\ - -^L 2 $ = (2.4) 

where L 2 is the usual squared angular momentum operator with eigenvalues L 2 = h 2 l(l + l). 

Since the problem is a spherically symmetric one, we put $ = exp((i/h)S(r, t))yj m (^, <p) 
where S is a function which will be expanded in powers of h. Substituting in the above 
equation, we obtain, 

[-(1 - 2M/rY l (d t S - eA t f + (1 - 2M/r) (d r S) 2 + (fi 2 + L 2 /r 2 )] 



+ - 



h 



-(1 - 2M/r)~ l d 2 t S + -^d r (r 2 (l - 2M/r))d T S + (1 - 2M/r)d 2 S 



(2.5) 



Expanding <S(r, t) = So(r, t) + (h/i) «Si(r, t) + (h/t) 2 S%(r, t) + ■ ■ •, substituting into Eq. (p.5|) 
and neglecting terms of order K/i and higher, we find to lowest order, 

- (1 - 2M/r)~ 1 (d t S - eA t f + (1 - 2M/r) (d r S ) 2 + (/i 2 + L 2 /r 2 ) = 0. (2.6) 

This is just the Hamilton- Jacobi equation satisfied by a charged particle of mass \i and 



charge e moving in the spacetime ( |2.1| ) and interacting with the potential A a as given in 
Eq. ( |2.3D . The solution to the above equation is 

So(r, t) =-et± [ dr(l - 2M/r)- 1 ^(e + eA t ) 2 - (1 - 2M/r)(/j 2 + L 2 /r 2 ) } (2.7) 



where e is a constant identified with the energy. 

The semiclassical kernel /C(r 2 , t%\ ri, ti) for the particle to propagate from (ti, ri) to (^2, ^2) 
in the saddle point approximation can be written down immediately: 
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IC(r 2 , t 2 ; n, t x ) = Afexp (~<%2, t 2 ; n, . (2.8) 

TV is a suitable normalization constant, and S is the action functional satisfying the classical 
Hamilton- Jacobi equation, namely 

S(r 2 , t 2 ; n, ti) = 5 (r 2 , i 2 ) - S (ri, h). (2.9) 

The sign ambiguity in So is related to the "outgoing" (d r S > 0) or "ingoing" (d r S < 0) 
nature of the particle. As long as points 1 and 2, between which the transition amplitude 
is calculated, are on the same side of the horizon (i.e. either exterior or interior to the 
horizon), the integral in the action is well defined and real. But if the points are located on 
opposite sides of the horizon, then the integral does not exist due to the divergence of the 
integrand at r = 2M. Therefore, in order to obtain the probability amplitude for crossing 
the horizon we have to give a prescription for evaluating the integral. The prescription used 
by Srinivasan and Padmanabhan 0, and adopted here, is to take the contour defining the 
integral to be an infinitesimal semicircle above the pole at r = 2M for outgoing particles 
on the left of the horizon and ingoing particles on the right. Similarly, for ingoing particles 
on the left and outgoing particles on the right of the horizon (which corresponds to a time 
reversed situation of the previous cases), the contour should be an infinitesimal semicircle 
below the pole at r = 2M. 

Consider, therefore, an outgoing particle (d r S > 0) at r = r 1 < 2M. The squared 
modulus of the amplitude for this particle to cross the horizon gives the probability to find 
the particle in the exterior of the horizon. The contribution to S in the ranges (r 1; 2M — 5) 
and (2M + 5, r 2 ) is real. Therefore, choosing the contour to lie in the upper complex plane, 

Sr ^ r * = -f™oJ 2M _ s T-2MTr dr+ (realpart) 

= i2-kM(s + A t (2M)) + (real part) (2.10) 

where A t is evaluated at r = 2M. The minus sign in front of the integral corresponds to 
the initial condition that d r S > at r = r\ < 2M . The same result is obtained when an 
ingoing particle (d r S < 0) is considered at r = r\ < 2M. The contour for this case must be 
chosen to lie in the lower complex plane. The amplitude for this particle to cross the horizon 
is the same as that of the outgoing particle due to the time reversal symmetry obeyed by 
the system. 

Consider next, an ingoing particle (d r S < 0) at r = r 2 > 2M. The squared modulus of 
the amplitude for this particle to cross the horizon gives the probability to find the particle 
in the interior of the horizon. Choosing the contour to lie in the upper complex plane, we 

get, 

Sr ^ r2 = -^L + s Y-2MTr dr+ (realpart) 

= -i2tiM(e + eA t (2M)) + (real part) (2.11) 

where as before A t is evaluated at r = 2M. The same result is obtained when an outgoing 
particle (d r S > 0) is considered at r = r 2 > 2M. The contour for this case should be in the 
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lower complex plane and the amplitude for this particle to cross the horizon is the same as 
that of the ingoing particle due to time reversal invariance. 

Taking the modulus square to obtain the probability, and substituting for A t from 
Eq. (p|), we get, 



probability for a particle\ I probability for a particle 

with energy e = exp (— 8tcMe + rre 2 ) I with energy e ] . (2.12) 

to be found at r > 2M J \ to be found at r < 2M 

This suggest that it is more likely for a particular region to gain particles than lose them. 
Therefore we must interpret the above result as saying that the probability of emission of 
particles is not the same as the probability of absorption of particles. Furthermore, this 
result implies that the flux of particle at infinity in this case would be greater by a factor 
exp(7re 2 ) than the flux of neutral particles of the same mass and spin. However, it should be 
noted that the increase factor is in fact negligiblely small when one considers the emission of 
electrons. In that case the numerical factor is equal to exp(7re 2 / (he)) ~ exp(7r/137) ~ 1.023. 
This result obviously cast doubt on the astrophysical importance of the electrostatic self- 
interaction. 

Note that the method presented here for the calculation of the modulation factor of 
Hawking's radiance is not specific to the electrostatic self-potential. It would also apply for 
an external electrostatic potential (as opposed to the internal electrostatic self-potential). 
Provided that the external potential is analytic in the neighborhood of the horizon, the 
thermal radiation is modulated by varying the magnitude of the external potential. The 
modulation factor is of the form exp(eA^ xt (2M)/Tfj), where the external potential, Af* is 
evaluated at the horizon and T# is the Hawking temperature. 

Recently van Putten [|TTJ discussed the magnetic analog of this phenomenon. He showed 



that a rotating black-hole produces electron-positron outflow when immersed in a magnetic 
field. The outflow is driven primarily by a coupling of the spin of the black hole to the 
fermionic wave-function. The source for the magnetic field is assumed to be exterior. He 
computed, in the WKB approximation, the superradiant amplification of scalar waves con- 
fined to a thin equatorial wedge around a Kerr black-hole and found it to be higher than for 



radiation incident over all angles. However, Aguirre 0] has recently presented calculations 
of both spin-0 (scalar) superradiance (integrating the radial equation rather than using the 
WKB method) and spin-1 (electromagnetic/magnetosonic) superradiance in van Putten's 
wedge geometry, and showed that in contrast to the scalar case, spin-1 superradiance is 
weaker in the wedge geometry. So that, as with the electrostatic self- interaction case, the 
astrophysical significance of the effect is questionable. 



B. The greybody factor 

It follows from the previous subsection that Hawking emission of charged self-interacting 
particles from the black-hole gives rise to a flux at infinity 

d 2 n 1 T . 

(2.13) 



dudt 2vr e 87rM (^- e2 /8M) _ 1 > 
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where uj is the energy of the particle at infinity, e is the charge of the particles, 8ttM is 
the inverse of Hawking temperature, and T is the relevant absorption factor, the so called 
"greybody factor". From the equation above it is apparent that the term — e 2 /8M serves 
as a chemical potential in the problem; intriguingly, this chemical potential is the same for 
both particles and antiparticles (usually chemical potential has opposite sign for particles 
and antiparticles). It may be seen that this chemical potential enhances the emission of 
(positive or negative) charged particles over that of neutral particles with the same mass 
and spin. Positive and negative charged particles are treated equally. 

It turns out that two different situations occur for uj > e 2 /8M and uj < e 2 /8M. To see 
this, expand the right hand side of Eq. fl2.13| ) around u = e 2 /8M. Then, we get that the 
flux at infinity is 

d 2 n 1 T 

(2.14) 



dudt 2vr 8irM(uj - e 2 /8M) ' 

It is now evident that for u > e 2 /8M the flux is positive. However for uj < e 2 /8M the flux 
seems to be negative, with a singularity at uj = e 2 /8M ! None the less, one shouldn't be 
perplexed by this. It is just that we have failed to take into account the energy dependence 
of the absorption factor T. Obviously, we must show that the the point of transition occurs 
at a zero of T, namely the expansion of the absorption factor around uj = e 2 /8M must begin 
with 

T « a{u, n, e 2 /8M)(w - e 2 /8M) + ■ • • , (2.15) 

where a(uj, fi,e 2 /8M) is finite function of its variables. 

Therefore, we turn to calculate T to first approximation. Making the ansatz $ = 
e~ lU}t Y l m (9 , ip)f(r)/r and substituting into Eq. (|2.4f) we obtain an effective Schrodinger equa- 
tion 

-h 2 ^ + V(r*)f = 0, 



V(r*) =(l-^)(/i 2 + ^+ / ' J2W 
2 



— \ L0 s-j 



e 2 M 
2r 



(2.16) 



Here r* = r + 2Mlog(r — 2M) is Wheeler's "tortoise" coordinate |I3"1| . We are interested 
in asymptotic solutions of the equation in the far region r* — > 00 (r — >• 00) and in the near 
region, r* — > — 00 (r — ► 2M). 

Taking the limit of V(r*) in the far region, r* — > 00 (r — > 00), we find that the equation 
has the form 



-^ 2 + (^ 2 -^)/ = O. (2.17) 



For n > uj this equation has exponentially decaying and diverging solutions. Obviously, 
the diverging solution is unacceptable on physical grounds, while the decaying solution is 
associated with a particle trapped in the effective potential well. For /x < uj Eq. ( [2.171) has 
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an ingoing and an outgoing wave solutions. Since we are actually dealing with a scattering 
problem, we concentrate on those solutions: 



f(r oo) = e-V^-^'-VR + ^t^-^/h ^ lg ) 

Here .A is a constant to be determined later. 

In the near region, r* —>■ oo (r — > 2M), we find that Eq. ( |2.16| ) has the limiting form 



72 -f 

h 2 -^-(uj-e 2 /8M) 2 f = 0. (2.19) 
dr* 



The Matzner boundary condition [14| that the physical solution be an ingoing wave, as 
appropriate to the absorbing character of the horizon, selects the solution of the above 
equation to be 

f(r -> 2M) = Be-t^/siiOrVR (2.20) 

where £> is a constant. 

A. and £> may be determined by matching / and /' of the solutions in the far and near 
regions at some point in the intermediate region 2M <C r <C oo. Doing so we find 

A = ^Hl^ e -^c 

uj - e 2 /8M + ^uj 2 - /i 2 
B = 2 v /u 2 -fi 2 e _|(^Z^- (a; _ e 2 /8M) ) c 

- e 2 /8M + v/cj 2 - /i 2 

(2.21) 

where r* m is the matching point. The absorption cross-section may now be obtained using 
the method of fluxes. The flux in our one-dimensional effective problem is 

T=^fd r f-c.c). (2.22) 

The absorption probability is the ratio of the incoming flux at the horizon to the incoming 
flux at infinity, 



^incoming 

4v47 ^ -(u-e 2 /8M). 



(u - e 2 /8M - ^uo 2 - fi 2 ) 2 

(2.23) 

The matching point r* m has disappeared from the expression for T. This means that at this 
order of approximation the matching point may be chosen arbitrarily as long as it is in the 
range 2M < r* m < oo. 

Note that for e 2 /8M w to this expression agree with our supposition (|2.15| ). The drastic 
difference between the two regimes, uj < e 2 /8M and e 2 /8M < uj shows up clearly if we 
consider the limit where the effective temperature of the black-hole tends to zero. If the 
self-energy of the particles is neglected, the rate of particles creation goes then to zero. In 
our system the rate goes also to zero if e 2 /8M < uj. But in the range \i < uj < e 2 /8M the 
rate of particle creation tends to —T/2tt. We are in the domain of black-hole superradiance 
of charged particles. 
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III. PAIR PRODUCTION NEAR A SCHWARZSCHILD BLACK-HOLE 



While Hawking radiance relates closely to dynamical spacetimes with horizons, pair cre- 
ation does not. One may inquire for the circumstances in which a neutral black-hole is 
involved in the spontaneous production of pairs of oppositely charged particles. Technically, 
pair production can take place when the conditions for the existence of a "generalized er- 
gosphere" (a region where orbits with negative total energy exist) are fulfilled. In this case, 
on the level of classical particles, a Penrose process [|T^| can take place. On the level of 



waves mechanics similar conditions must be fulfilled in order for "superradiant" scattering 
of waves obeying the Klein-Gordon equation to occur. The corresponding phenomenon that 
occurs then is the so called Klein paradox . However, all this is known not to happen 
when the black-hole is a neutral static and isolated one. But, it turns out that once the 
black hole is allowed to interact with an exterior charged test object, this may give rise to 
pair production. Technically, this means that the conditions for the occurrence of a Klein 
paradox are obeyed. Then it is possible to calculate in a W.K.B. approximation the trans- 
mission coefficient through the Klein region separating the positive from the negative states 
of the corresponding Klein-Gordon equation. This gives the probability for pair creation. 

Using the Hamilton- Jacobi equation Eq. ( |2.6|) , one may derive the equation of motion of 
a classical particle of mass \i and charge e on the background metric (|2.1|): 



dr\ 2 
dr) 



(e-S + (r))(e-S4r)), 



Here r and r are the radial coordinate and proper time of the particle respectively, and S + 
and £_ are the effective potentials for the positive and negative energy solutions, respectively. 
The classical bound states in the potentials £± are the classical limit of the "resonances" of 
a quantum field satisfying the Klein-Gordon equation (|2.2| ) written in the given background 
metric. Positive energy states, e + > £+, correspond to a positive probability density, and 
therefore describe particles of energy e. Negative energy states, e_ < £_ correspond to a 
negative probability density and therefore describe antiparticles of energy —e. When there 
is a crossing of the e + and the £_ states, the probability density has a variable sign. We are 
then in a Klein paradox region. 

The transmission coefficient T 2 through the potential barrier separating the positive 
from the negative energy states is proportional to the probability for an incident particle to 
create a pair of particles. The transmission coefficient can be computed using the W.K.B 
approximation to Eq. (|2.16|) : 

T 2 = exp(-C), 



C = 2 I ^mt;^ (3.2) 



«2 



Here ( may be identified as the opacity of the barrier against pair creation, cti and ai are 
the zeroes of the effective potential V(r) defined in Eq. ( [2.16| ). 



S 



What then is the energy range for which a Klein region exists? The Klein region is 
defined by 



min£ + < e < maxS. 



(3.3) 



Now, £_ is a monotonic decreasing function of r. It attains its maximum on the horizon: 
max£_ = £_(r = 2M) = e 2 /8M. Its asymptotic value at infinity is —[/,. £ + has a positively 
diverging derivative on the horizon where it attains the same value as and an asymptotic 
value at infinity, /i, which it approaches with a positive slope. This obviously implies that 
somewhere in the intermediate region between the horizon and infinity, there is a point 
where S + attains a minimum, e.g. a potential well. The location of this minimum may be 
found by calculating the roots of the equation d r £ + = 0. There are two physical roots. One 
determines the location of the potential minimum while the other determines the location of 
the potential maximum. Now, the effect of angular momentum may be ignored, since this 
will only strengthen the potential barrier leading to a smaller particle creation rate. Doing 
so, we may ask what are the allowed values of 7 = e 2 /8M/x for the existence of a Klein 
region. Obviously, for 7 = (e = 0) which corresponds to the usual neutral Schwarzschild 
spacetime, there is no level crossing, no Klein region, and thus no particle production. This 
indicates that we should actually look at the other extreme, where 7 is large. However, for 
an electron e 2 /(G/z) ~ 10 15 g. This means that for 7 to be large we must look at black-holes 
of mass not much bigger than 10 15 g: the arena of mini black-holes. Keeping this in mind, 
we find the minimum of the effective potential £_ for large 7 which finally gives the energy 
range for the existence of a Klein region: 



The function in the square brackets in the leftmost side of the inequality above tends rapidly 
to 1 as 7 is increased (for 7 = 1 it is equal to 0.926; for 7 = 5 it is equal to 0.987). Therefore, 
we conclude that in practice, the energy range for the existence of a Klein region is 



This energy range is the same as the one in which the black-hole starts to superradiate (see 
the end of the previous section). 

It should be emphasized that unlike the case where the background carries a definite sign 
of electric charge (like in the Kerr-Newmann spacetime), here the background is neutral. 
Therefore, it would seem that the black-hole should absorb statistically equal amounts of 
particles and antiparticles. Hence, it should remain neutral. In reality, a screening effect 
should take place. Since the point charge which creates the electric field should repel particles 
(or antiparticles) with the same sign of charge as its own, and attract their counterparts, a 
cloud of charge would form around it. This in turn should screen the point charge and further 
lower the pair creation rate. Furthermore, this charge segregation may alter the probabilities 
for assimilation of particles (antiparticles) by the black-hole, so that assimilation of particles 
with sign like the sign of the test charge would be more probable. Thus, the black-hole may 
become charged after all ! 




(3.4) 
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Unfortunately, calculating the corresponding Debye length of the problem is a formidable 
task and falls outside the scope of this work. Surely, the cloud of charge may be considered 
as a neutral gas consisting of charged particles which Coulomb-interact. Then, in order to 
solve the problem one must first solve Maxwell equations on the curved background with a 
source term representing the distribution of charge in the cloud. In the case of a completely 
ionized gas or plasma, this source term is given by a sum of Boltzmann's factors, one for 
each kind of particles, each multiplied by the charge carried by each kind of particles ]TT 



However, this is valid only if the pairs may be considered to be thermally distributed, which 
is not generally true in the case at hand. 



IV. WHERE DOES THE ENERGY COME FROM? - A SPECULATION. 

In the previous section it was shown that the system particle-black-hole loses energy in 
the process of pair creation. A simple question may be raised - what is the source of energy 
carried by the pairs? The corresponding problem of the possible sources of electromagnetic 
energy radiated away by an accelerated charge in flat spacetime, troubled and still troubles 
scientists [|18|] (it was named as the "Energy Balance Paradox"), and several answers were 
suggested. Leibovitz and Peres ITS] suggested that there exists a charged plane, whose 



charge is equal and opposite in sign to the accelerated charge, and that it recesses with the 
speed of light in a direction opposite to the direction of the acceleration. The interaction 
between this charged plane and the accelerated charge supplies the energy carried by the 
radiation. Another suggestion, by Fulton and Rohrlich []20|], is that the energy radiated is 



supplied from the self-energy of the charge. In the problem considered here there is a third 
energy source-the black-hole. Here, we show that on thermodynamical grounds the option 
that the black-hole lose mass during the process is not possible for massive black-holes, and 
discuss the proposition that the self-energy of the test charge is the energy source. 



A. Thermodynamical arguments 

First assume that the black-hole loses energy given by —AM. This can be approximated 
by —AM m 2eN, where N is the number of pairs, and e is some mean energy carried away 
by the pairs. In the nonrelativistic regime, e « // where 2/i is the pair rest-mass. In the 
relativistic regime the rest-mass of the particle can be neglected, so /x <C e. Now, since the 
black- hole is static and neutral, its entropy is given simply by Sbh — 4itM 2 . Therefore, 
during the process the black-hole changes its entropy by ASbh = SttMAM. Combining 
the two results, we find that 

AS BH = -!6nMeN. (4.1) 

Now we make use of the generalized second law of thermodynamics (GSL). For the GSL to 
hold, the entropy carried out by the pairs must at least compensate for the entropy lost by 
the black-hole, namely 

< AS wor id = ASbh + *AS pa i rs . (4-2) 
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Now, if the created particles are considered to be nonrelativistic, the entropy they carry is 
never far from the number of particles involved. Thus, 



ASpair, ~ VN, (4.3) 

where rj is a proportionality constant of order unity. The same is known to be true in 
the other extreme when the created pairs are assumed to be relativistic. For example, for 
black body radiation 77 = 27r 4 /(45C(3)) ~ 3.6, where ((z) is the Riemann zeta function pT 



Similarly, if the duration of the pair production process is long, so that the pairs are allowed 
to thermalize, they may be considered as particles obeying Fermi-statistics with vanishing 
chemical potential [|17 | . The specific entropy is then 77 = S/N = 147r 4 /(135C(3)) pa 8.4. 
Substituting ASbh and AS pairs into the GSL we find 

- 16nsMN + r)N>0. (4.4) 

Hence 

fiM < eM < — = 0(1). (4.5) 
16n 

The conclusion must be that in the 1 <C fiM regime, the black-hole cannot be the dominant 
energy source; the black-hole is just too cold! Accordingly, the energy must come from 
somewhere else. A further conclusion is that the black-hole may not lose entropy during the 
process: the black-hole is involved in an adiabatic process making its horizon area invariant 

Actually, in the last derivation we have implicitly assumed that the pairs are produced 
at a large distance from the horizon. It turns out that the GSL even strongly forbids pair 
production by the black-hole taking place at close proximity to the horizon. There the 
energy of the pairs as measured locally is dominated by the electrostatic self-energy, and in 
fact diverges. To see this we note that if the constituents of a particle-antiparticle pair are 
considered to be quasistatic, then their conserved energy as measured at infinity, e, amounts 
to energy invested in rest-mass plus electrostatic self-energy: 

e = (1 - 2M/r)e locol 
= (1 - 2M/r) f/i + —^L—\ , (4.6) 



2r(r - 2M) 

where r is the Schwarzschild coordinate of the particle (see Eq.( |4.12j ) below). The first strik- 
ing thing apparent from the expression for Si oca i is that it diverges as the particle approach 
the horizon, r — > 2M. Therefore, if the pairs are located in the close proximity of the horizon 
then the dominant part of their energy lies in electrostatic self-energy. In the other extreme, 
when the particles are located far away from the horizon, their electrostatic self-energy is 
small compared with the rest-mass energy, until it vanishes at infinity. The point of transi- 
tion, f, from electrostatic self-energy dominated region to rest-mass dominated region is set 
by the condition 

_ q 2 M 
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Thus f = 2M/(e 2 /8M/i), where f(x) = ^x + 1/4 + 1/2. Note that /(0) = 1 - no elec- 
trostatic self-energy. As was assumed in the previous section, the pair production rate is 
exponentially small. Therefore in order to produce a non-negligible number of pairs, the 
exponent must be of order unity, e.g. we are limited by the condition 1 < e 2 /M/z (see 
Eq. (|3.5|) ). Taking this into consideration we note that f(l) ~ 1.62. Hence, in our approx- 
imation, the electrostatic self-energy is the dominant part of the particles energy only in a 
narrow region around the horizon of width, r — 2M ~ 0.62 x 2M. Pair production in that 
region is highly improbable on account of the GSL. 



B. Dynamical arguments 

Realizing that the black- hole may not be the major energy source, we turn now to look at 
the dynamics of the test charge as a second candidate. We begin by considering the motion 
of a test particle of mass m and charge e. Its motion, were it subject only to gravitation and 
electromagnetic influences, would be governed by the Lagrangian 

L = -m J yj — g a p u a u@ dr + e J A a u a dr, (4.8) 

where x a (r) denotes the particle's trajectory, r the proper time, and u a = x a = dx a /dr, and 
A a means the background electromagnetic 4-potential evaluated at the particle's spacetime 
position. Recalling that g a pu a VtP = —1, it follows from the Lagrangian that the canonical 
momenta are p a = 5C/5u a = mgapu 13 + eA a . The stationarity of the envisaged background 
means there is a timelike Killing vector £ a = {1,0,0,0}. The quantity 

e = -p a C = -mg tl3 u 13 - eA t , (4.9) 

corresponds to the usual notion of energy as measured at infinity. Its first term expands 
to m + ^m(dx/dt) 2 in the Newtonian limit. The second term, —eA t , is thus the electric 
potential energy. 

Varying the Lagrangian Eq. ([4.8|) with respect to u 13 gives the equation of motion of the 
particle 

m=--— = eF.V 3 , (4.10) 
dr 

where F a p = Ap. a — A a -^. The proper time derivative of e, may be calculated as follows: 

t = - -e (m^ + e2£) - eAa Bl. (4.11) 



dr \ dr dr J dr 

where we have used the fact that the proper time derivative of the timelike Killing vector, £ Q , 
along the trajectory of the particle is always perpendicular to the trajectory, {u a D^ a / dr = 
0). 

Now, If the test charge is assumed to be quasistatic (supported by some mechanical 
apparatus), then its 4- velocity is given by u 13 w ((-gu)- 1 ' 2 , 0,0,0). Substituting for A a 
from Eq. ( |2.3| ) into Eq. ( |4.9| ) we find 



12 



m(l -2M/r) + e 2 M/(2r 2 ). (4.12) 



Making use of the equation of motion ( |4.10| ), it is easy to show that e vanishes. Even 



if we go beyond the quasistatic approximation and assume that the particle has a small 
but non-vanishing radial velocity, e is still conserved. Does it mean that the system cannot 
radiate? We intend to show now how the equation of motion should be modified to account 
for the irreversible processes involved in pair production. 

First, we define a 4-momentum rate of radiation by lZu a . In the case where the radia- 
tion is electromagnetic TZ is given by the relativistic generalization of the Larmor formula 
TZ = (2/3)e 2 a a a a ; a a is the acceleration. In the problem considered here, a simple relation 
between TZ and the dynamics of the test particle is unknown. However, one may approximate 
the rate of energy loss due to pair production using methods to be discussed below. 

TZu a , being a loss, should be subtracted from the right hand side of the equation of motion 
( |4.10| ). One might expect that in this way we have correctly accounted for the momentum 
and energy loss due to pair production. Unfortunately, this equation is inconsistent with TZ 
being positive definite; multiplication of the modified equation of motion by u a , and using 
the normalization condition, u a u a = — 1 yields TZ = 0! A term must be missing. To supply 
it we write 

m—— = eFZu 13 - T a , 

GST P 

Y a = TZu a + S a . (4.13) 

r a here may be understood as a 'frictional force', and S a is to be specified below. 

There remains the question of what is the origin of that friction? One possible answer 
is to view the test particle as a Brownian particle interacting with a quantum field assumed 
to be in the vacuum state f24|. At zero temperature, even though the thermal fluctuations 



are absent, the quantum field still possess vacuum fluctuations. Obviously, the particle 
cannot keep accruing energy from the fluctuations present in the surrounding environment. 
Therefore, there should exist a mechanism for the particle to dissipate its energy so that 
it reaches equilibrium with the environment. Treating the quantum field as a classical 
stochastic variable it is possible to take the approach of Langevin who suggested, early this 
century, that the force exerted on the particle by the surrounding medium can effectively 
be written as a 'rapidly fluctuating' part and an 'averaged out' part which represents a 
frictional force experienced by the particle. The presence of the frictional force implies the 
existence of processes whereby the energy associated with the particle is dissipated to the 
degrees of freedom corresponding to the surrounding medium. 
Multiplying Eq. fl£I3D by u a now yields 

T a u a = TZ = S a u a . (4.14) 

We have obtained a constraint over S a . It is reasonable to assume that S a is a function of 
the velocity and its derivatives. Given that the variation of the velocity is small, we may 
expand 

S Q = C u a + da a + C 2 a a + C 3 d a • • • , (4.15) 

where {Cj} are proportionality constants with corresponding dimensions of energy x time 1 ' 1 . 
The first term has to vanish since it is already included in lZu a . The second term may be 
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accounted for by mass renormalization by the rule, m — > m + C\ (see the equation of motion 
( |4.13| )). It turns out that the proper time derivative of the acceleration, a a , is the lowest 
derivative of the velocity allowed. Accordingly we set S a = C2 d a . 

In conjunction with the constraint (]4 . 1 4j) , we must go beyond the static approximation. 
For if the particle is static then uqu° = —1, and 1Z = u°So- Thus , a straightforward 
calculation of e using the modified equation of motion (|4.13 ), gives 



e = -£ a r a = -(u u°S + So) = 0, (4.16) 

regardless of the properties of S a . Accordingly, we set 

u a = ((-9tt)- 1/2 (1 + g r M 2 /2) , ±8u, 0, 0) , (4.17) 

where the velocity correction, Su, is assumed to be time independent and small. Calculating 
a a to the lowest order in Su and its derivatives, and substituting into the constraint Q4.14|) , 



we obtain a first order, non linear, ordinary differential equation for Su(r), whose solution 
is: 



Su{r) = ± 



MCo 



-9ttf 



1/4 / r r \ 1/4 



r 2 lZdr . (4.18) 



Su scales with (1/M) 1//4 . Hence, for massive black-holes Su is small. Moreover, Su is pro- 
portional to the 1/4-th power of the energy dissipated along the trajectory of the particle, 
hence it depends on the history of the particle. TZ must drop at least as fast as 1/r 4 for 5u 
to converge as r — > 00. In fact, it was already assumed in Sec.|nT| that 1Z is proportional 
to the transmission probability through the Klein region which is exponentially small. This 
however determines only the energy scale over which pair production may take place. It 
does not set the functional dependence of 1Z at the position of the test particle. 

To resolve this we turn to Schwinger's approach [p5|, who showed that the probability 



for pair creation per unit time per unit volume by a constant electric field is 

qE \ 2 *^ 1 fc^M 2 . 

7 fc=i 

Multiplying this with some mean energy carried by the pairs gives an approximation for 
1Z. Although (|4.19| ) was originally derived by calculating the imaginary part of the effective 



Lagrangian for the Dirac field of rest-mass /i and charge q in a prescribed constant elec- 
trostatic field E in flat space time, we adapt this result to the corresponding problem in 
curved spacetime by substituting local expressions for the energy and electric field. We use 
the formula for the repulsive self-force as measured by an instantaneously comoving, freely 
falling observer, at the position of the test particle 0], to define the effective electric field E: 

e 2 M 

E self = eE = —. (4.20) 

This repulsive force is peculiar to charged test particles. Since the black-hole is uncharged, 
this must be interpreted as arising from the test particle's electrostatic self-interaction. It 
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vanishes as M vanishes, indicating that the effect is induced by the black-hole's spacetime 
curvature. 

Substituting for the effective electrostatic field in the formula for the rate of pair pro- 
duction, Eq. (|4.19p ) we find that 71, vanishes exponentially fast as r — > oo. Finally, putting 
everything together we find that 

£ = - = 0{n*'% (4.21) 

hence, the change in the energy of the test particle is negligiblly small regardless of the sign 
of the velocity. 

To summarize, as the particle is slowly lowered towards the black-hole (or pulled back) by 
the mechanical apparatus, additional work must be done against the frictional force induced 
by vacuum fluctuations. The extra energy invested in moving the particle is then dissipated 
away as pairs of massive charged particles. 



V. SUMMARY AND ASSESSMENT 

It was demonstrated how Hawking's radiance form an isolated neutral black-hole is mod- 
ified on account of the electrostatic self-interaction of charged particles. Once the temper- 
ature of the black-hole is high enough so that the lightest massive charged particles are 
emitted, the thermal radiation of charged particles emitted by the black-hole is increased 
with respect to the thermal radiation of neutral particles with the same mass and spin. This 
is a direct consequence of the inclusion of the self-interaction into the analysis. 

An interesting consequence of this conclusion is that an external electrostatic potential 
(as opposed to the internal electrostatic self-potential) can also be used to modulate Hawk- 
ing radiance. Provided that the external potential is analytic in the neighborhood of the 
horizon, the thermal radiation is modulated by varying the magnitude of the potential. The 
modulation factor has the form exp(eAf xt (2M) /Th), where the external potential, Af^ is 
evaluated at the horizon and T# is the Hawking temperature (see Eq. (|2.12|) ). 

The possibility to modulate the thermal emission from the black-hole has some very 
interesting consequences. First, assume that one applies an external repulsive electrostatic 
field, opposing the gravitational pull of the black hole. Then the probability to propagate 
from the interior to the exterior of the horizon for a charged particle with energy below 
eA^ xt (2M), would be greater than the probability for the inverse process. If now the applied 
electrostatic field is attractive (acting in the direction of the gravitational pull), then it 
would serves as a high pass filter, suppressing the emission of charged particles with energy 
below eAf xt (2M). This enable us to control the energy range over which the black-hole 
superradiate charged particles! The similarity between this effect and the phenomenon of 
thermionic emission is manifest. 

Another issue that deserve further investigation is the calculation of the rate of pair 
production. It was assumed in Sec.[TTT] that the transmission coefficient, T 2 , through the 
potential barrier separating the positive from the negative energy states is proportional to 
the probability for an incident particle to create a pair of particles. The Klein region for 
the problem (including the effect of self-interaction) was found and shown to correspond 
to the energy range for black-hole superradiance of charged particles. Thought T 2 may be 
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calculated numerically, it would be profitable if one could find an analytic approximation for 
T 2 , and show that the result converges to the result obtained using Schwinger's approach 
26| . The problem is that, Schwinger's approach was originally formulated in flat spacetime, 



and the formulation of this approach in curved spacetime may very much prove to be an 
uphill task. So that a way to bypass these difficulties is much needed. 

Finally, the problem of energy source for the pair production process was discussed. It 
was shown, that on thermodynamical grounds, it is not possible for a massive black-hole 
to lose mass during the process. This is just to say that the black- hole is too cold. More 
precisely, the entropy outflow from the system is too low for the generalized second law to 
hold. We thus turned to explore the dynamics of the test charge as a second candidate. 
It was speculated that vacuum fluctuations of a quantum field interacting with the test 
particle may be involved in the process of pair production. These vacuum fluctuation, 
induce random motion that the particle undergoes, and an averaged-out force that enters 
into the equation of motion as a friction term. This friction term is a manifestation of the 
dissipation mechanism by which energy is given off in the form of massively charged particles. 
That being the case, a relation between the friction term and the rate of energy dissipation, 
was found. As could be anticipated, a static system can not radiate. Accordingly, going 
beyond the static approximation, it was assumed that the test particle has a small (but 
non-negligible) radial velocity. Then, the functional dependence of the particle velocity on 
the rate of energy dissipation was determined. It was shown that the velocity scales inversely 
with the black-hole mass and proportional to the 1/4-th power of the energy dissipated along 
the trajectory of the particle, hence it depends on the history of the particle. 

The picture that seems to arise is that as the particle is lowered towards the black hole, 
or pulled away, the mechanical apparatus supporting the particle is doing work in changing 
the particle's energy to the value appropriate to the new location. However, as the particle 
moves, it interacts with the vacuum fluctuations in the medium which have the effect of 
inducing frictional forces. If the particle is assumed to move in a constant velocity, then 
these frictional forces must be overcome by the investment of additional work on part of 
the mechanical support. The extra energy is then dissipated away in a process of pair 
production. Obviously for the establishment of this picture, further study of the relationship 
between vacuum fluctuations in curved spacetime and friction is in order. 
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